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We study the non-ideal effects arising due to viscosity (both bulk and shear), equation of state 
(e 7^ 3P) and cavitation on thermal dilepton production from QGP at RHIC energies. We calculate 
the first order corrections to the dilepton production rates due to shear and bulk viscosities. Ignoring 
the cavitation can lead to a wrong estimation of dilepton spectra. We show that the shear viscosity 
can enhance the thermal dilepton spectra whereas the bulk viscosity can suppress it. We present 
the combined effect of bulk and shear viscosities on the dilepton spectra. 
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I. INTRODUCTION 

QGP formed in the relativistic heavy ion coUider (RHIC) experiments is considered as the most perfect fluid in 
the nature [IH3]. Experiments point towards a very low value for the shear- viscosity to entropy density rj/s ^ I/Att- 
the KSS limit [1], for this strongly coupled matter formed at RHIC [5]. Generally the effect of bulk- viscosity C is 
neglected for a system obeying a relativistic equation of state (EoS). This is because ( scales like — |, where is 
the speed of sound. Thus for a system with ideal EoS (e « 3P), where e is the energy density and P is the pressure 
= 1/3 = c2), the effect C is neglected [6] . However in the temperature regime close to the the critical temperature 
Tc, the bulk viscosity may not be a negligible quantity. In fact, the recent lattice QCD results indicate that the 
quark-gluon matter EoS departures from the ideal EoS description near critical temperature Tel?]. The ratio of the 
bulk viscosity to the entropy density s calculated from the lattice data indicate C/s too has a peak and >> 77 at 
the regime T k, Tc [8l-fT3]. Second order relativistic hydrodynamics models have been used to describe the role of 
viscosities dynamics on the evolution of the system |14H20| . Using one dimensional hydrodynamics role of bulk viscosity 
in heavy ion collisions is recently analvsed 21, ,22]. Incorporation of bulk viscosity into heavy ion scenario can bring 
in interesting phenomena of cavitation. Since bulk (shear) viscosity reduces the longitudinal pressure of the system, 
during the evolution of the system with sufficient values of the viscosities, effective pressure of the fluid can become 
zero causing cavitation. Cavitation leads the fluid to break apart into fragments, making further hydrodynamical 
description invalid In addition the viscous effects can modify the temperature profile and thereby it can 

change the particle distribution functions of the plasma(25]. Using kinetic theory methods one can include these 
corrections in the distribution functions and this may have observable consequences in the observables |26ff29] . 

Thermal photons and dileptons are among the most promising probes of the hot and dense matter created in 
relativistic heavy ion collisions [301 ISI]- As their mean free path is larger than the transverse size of the fireball, they 
can escape from the system and there by provide information about the thermodynamic state and space-time history of 
the matter created in heavy ion collisions [32 . Production rates of these probes (particles) depend on the temperature 
of the system and by knowing the appropriate initial conditions time evolution of the temperature of the system can 
be obtained by using the equations of the hydrodynamics. Once the temperature profile is obtained, the calculation 
of the thermal spectra can be done by evaluating the cross-section of the underlying scattering processes. We refer 
readers Refs.[33 37J for excellent reviews on the subject. Thermal photon production have been studied under various 
conditions by several authors [55ti43j . Thermal photons from quark-gluon plasma (QGP) in the presence of shear 
viscosity was studied recently in Refs.f44-'3S] and they were proposed as a tool to measure the shear viscosity of the 
matter formed in the heavy ion collisions [44, t45]. As argued earlier in the temperature regime T ^ Tc, the finite 
bulk viscosity may significantly influence the hydrodynamic evolution of the system [2TJ [22]. Recently the role of these 
non-ideal effects due to EoS, bulk viscosity and cavitation were considered by us in thermal photon production|47|. 
We showed that all these effects can alter the particle spectra in a significant manner. It was seen that if the effect 
of cavitation were not included properly one will end up with erroneous estimates of the particle production rates. 
Thermal dilepton production using the equations of ideal hydrodynamics is well studied by many authors [48, -.50] . 
However, only recently the thermal dilepton-production from QGP in the presence of shear viscosity was studied[5T]. 
In the present work we investigate the effects of finite bulk-viscosity on the thermal dilepton production from QGP. 
The main source of thermal dileptons is from the quark-anti-quark annihilations: — >■ 7* — > l^l^ . The cross-section 
of this lowest order process is well known [52]. There are other higher order processes which may also contribute in 
thermal dilepton production |531 [51] . however we are not considering them in this present analysis. Thermal dileptons 
from the annihilation process is dominant in the window of intermediate invariant mass and transverse momentum of 
the lepton pair 1 < M, < 3 GeV [55[[56] . 

In Ref.[5T] authors has studied the role of shear viscosity in thermal dilepton spectra. However we note that they 
had used an ideal EoS for the calculation and the effect of bulk viscosity was not considered. It has been shown in 
our previous work that bulk viscosity plays a dual role in heavy ion collisions, on the one hand it enhances the time 
by which system cools down to Tc and on the other hand it can make the hydrodynamic treatment invalid much 
before it reaches Tc. Non-ideal EoS also makes the system spend more time near Tc thus increasing the thermal 
particle production '47'. In this work we are studying the effect of non-ideal EoS, bulk viscosity and cavitation on 
the thermal dilepton production from the QGP. Firstly, by taking the viscous modified distribution functions using 
the 14-moment Grad's method results, we calculate the first order correction due to both bulk and shear viscosity in 
the dilepton production rate. We use a recent lattice QCD calculation result for non-ideal EoS. We take ij/s = 1/Att 
in our analysis and we use lattice QCD result for (/s following Ref.[22 . By using second order relativistic causal 
dissipative hydrodynamics we analyse the evolution dynamics of the system. Where we treat the expanding system as 
one dimensional boost invariant flow. Since boost-invariant hydrodynamics leads to an underestimation of the effects 
of bulk viscosity, our dilepton spectra acts as a conservative estimate of the effects. At early stages of the expansion, 
transverse flow can be neglected. Eventhough we are not including transverse flow, we believe its effects could remain 
small as cavitation can reduce the hydrodynamical evolution. 
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II. THERMAL DILEPTON PRODUCTION RATES IN QGP 

In QGP the dominant mechanism for the production of thermal dileptons comes from qq annihilation process 
qq ^* — !■ l^l^ . From kinetic theory rate of dilepton production (number of dileptons produced per unit volume 
per unit time) for this process is given by 



where pi,2 — (-©1,2, Pi 2) is the four momentum of quark or anti-quark with Ei^2 = yPi 2 + "^g — |Pi 2 1 neglecting 

the quark masses. Here — {Ei + £'2)^ — (Pi + P2)^ is the invariant mass of the virtual photon. The function 
f{E,T) = 1/(1 + e^^"^) is the quark (anti-quark) distribution function in thermal equilibrium and g is the degeneracy 

factor. Vrei = ^AElE^''^ ^ 2E1E2 rclativc velocity of the quark-anti-quark pair and cr{AP) is the thermal 

dilepton production cross section. The cross-section a^M'^) in the Born approximation is well known: g'^a{M^) — 

^'^'"^ ^I^i^^^^^" ^^'^ ^i^'^ ^/=2 and = 3, we have M'^g'^a{M'^) = ^a^^. Since we are interested in the rate for 
a given dilepton mass and momentum, we write 

d^p-J (2^)3 (2^)3 2E,E2 ^ ^P~P^-P'^ (2) 

where p — {po = Ei + E2, P = Pi + P2) is the four momentum of the dileptons. At the present case we are interested 
in the invariant masses larger compared to the temperature i.e.; M ^ T, in this limit we can replace Fermi-Dirac 
distribution with classical Maxwell-Boltzmann distribution, 

/(i?,T)^/o = e-^/^. (3) 



III. VISCOUS CORRECTIONS TO THE DILEPTON PRODUCTION RATES 

Viscosity effects modify the particle production spectra in two ways. Firstly it modifies the temperature and 
secondly through the corrections in distribution functions |47!. The first effect effect is incorporated when we calculate 
the temperature as a function of time by solving dissipative hydrodynamics. This will be done in the next section. 
Once we include bulk viscosity, a novel phenomena like cavitation can arise which can alter particle production spectra 
considerably J47^. In the present section we concentrate on the second effect. To calculate the viscous modifications 
to the distribution function as a function of the momentum, we need to use the techniques of relativistic kinetic 
theory [2^. Let us write the modified distribution function as f = fo + Sf, with viscous correction 5f = Sf,j + Sf(^, 
where Sfrj and 6f(^ represent change in the distribution function due to shear and bulk viscosity respectively. We can 
calculate these corrections using 14-moment Grad's method as done in Refs.|26| , [27 | l47]. Now modified distribution 
function from corrections due to shear (77) and bulk {() viscosity (up to quadratic order of momentum) are given 
bygT] 

f{p) = Mp)(i + ^py'Vi^Ufs) + ^^pVA^^e). (4) 

where s is the entropy density and the operators are defined as A"^ = g"^ — u" , — Aajjd'^ , Q = V^u" and 

In order to compute the effect of viscosity on the production rate, we substitute equation Q in dilepton rate 
equation Thus keeping terms up to the second order in 77/s and C/s, the dilepton production rates can be written 



dN diV(o) dAf('') d7V(^) 

(5) 



d'^xd^p d^xd-^p d-^xd^p d^xd'^p' 
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(6) 
(7) 
(8) 



The first term (given by equation ([6|)) is the one without any viscous corrections {ideal part) and is weh known 

d'^xd'^p 2 (27r)5 

The first order correction to the rate due to shear viscosity- given by equation ([?]), is calculated in Ref. [5T] and the 
final expression is 



(9) 



d^xd^p 2 (27r)5 ^ 3 



2r3 



(10) 



Let us next proceed to estimate the correction to the rate due to bulk viscosity from equation (|8|. We can write 



dN<-0 



d^p. 



-{Et + E2)/T 



d'^xd^p J (27r)6 
where we have represented 



5{p 



,-E,-E2) = m^I-^p)^ 



5 T3 



./^e, (11) 



(12) 



Now we write the second rank tensor J"*^ in the most general form constructed out of and p": 

= ao5"'^ + aiu^'u^ + a2pV + ci^{u°'p^ + u^p°') ■ (13) 

Note that because of the identity u^Aap = 0, the coefficients of /"^ which are going to survive after contraction 
with Aq,^ are oq and 02- We construct two projection operators to get these coefficients, i.e.; Q^jjl"^ = oo and 

Oa,3^"'^ =«2, so that 



The expressions for the projection operators in the local rest frame of the the medium {u" = (1, 0)) are 

QIb = WrJ^[\P\'^9aP + M'^UaUfS + PaPf3 - 2poUaP/3] , 



2|pr 

1 



Qa/3 = ^|^[|Pl 5a/3 + (3Po ^ |Pl )UaUp + SpaPp - GpoUal 



(14) 

(15) 
(16) 



With the help of definition of /"^ i.e.; equation (12 1, we can calculate {Qj^J'"') and {QIJ^""). We now write the 
final expression for the first order correction due to bulk viscosity in dilepton rate 



d'^xd'^p 2 (27r)5 ^ 



3 V 5 2r3 



5 4r3 



(17) 



where we have used the identity A|^ = 3. 

The total dilepton rate, including the first order viscous corrections due to both shear and bulk viscosity is obtained 
by adding equations (10 1 and (17). 

Apart from rates as function of four momentum of the dileptons we will be interested in particle production as 
a function of invariant mass (M), transverse momentum (pt) and rapidity (y) of the dilepton pair. This can be 
obtained from equation ^ by changing the variables appropriately (33) . 

dN 1 dN 



d^xdM'^d^prdy ~ 2 d^xd^p 
_l W 
" 239^ 



■Po/T 



(18) 
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IV. VISCOUS HYDRODYNAMICS AND CAVITATION 



In order to study the dilepton production from the QGP formed in heavy ion coUision we need to understand the 
evolution dynamics of the system. By obtaining the temporal distribution of temperature and information about 
viscosity coefficients we can study the dilepton spectra. 

We study the QGP formed in high energy nuclear collisions using causal dissipative second order hydrodynamics of 
Israel-Stewart [57] . with the expanding fireball treated as having one dimensional boost invariant expanding flow|58j. 
We use the parametrization of the coordinates t = t coshyy^ and z = t sinh?7<,, with the proper time r = y/t'^ — 
and space-time rapidity rjs — '"-[f^f ]• Now the four velocity can be written as 



(cosh 77s, 0, 0, sinh rjs). 



(19) 



With this second order theory (For more details on this theory and its application to relativistic heavy ion collisions 
we refer [ITl [321 [30] ) the equations dictating the longitudinal expansion of the medium are given by 



de 










f P + 






a$ 




2 1 




h 


3^ 


an 


n 


1 




m 





n- $) 
1 



4tv 
3r 



A $2' 
2r/2 



(20) 
(21) 
(22) 



The effects due to shear and bulk viscosity are represented via $ and 11 respectively and they alter the equilibrium 
pressure. The first equation is the equation of motion and the other two equations: (21 & 22) are evolution equations 
for $ and 11 governed by their relaxation times and Tn respectively, 
the relaxation time Tn = C /^o and = 
and Ai IS5HH7I: = and Ai 



27rT 



Apart from these three equations (20 



The coefficients /3o and fi2 are related with 
277/32- We use the N = A supersymmetric Yang-Mills theory expressions for 
- ^^and we take r^(T) = Tn(r) following Ref.[21 . 

22), we need to provide the EoS to study the hydrodynamical evolution 



of the system. We use the recent lattice QCD result of A. Bazavov et aZ.|7] for equilibrium equation of state (EoS) 
{non-ideal: e — 3P ^ 0), which becomes significantly important near the critical temperature. Parametrised form of 
their result for trace anomaly is given by 



e-3P 

y4 



= 1 



1 



1 -I- exp 



rp2 



(23) 



where values of the coefficients are ^2 = 0.24 GeV2, d^ = 0.0054 GeV^, ci = 0.2073 GeV, and ca 0.0172 GeV. The 
functional form of the pressure is given by |7] 



PjT) PjTo) 



dT 



„ £ - 3P 



(24) 



with Tq — 50 MeV and P{Tpi) = (22J. From equations (23 1 and (24) we get e and P in terms of T. A crossover 
from QGP to hadron gas around the temperature 200-180 MeV is predicted by this model. Throughout the analysis 
we keep the critical temperature to be 190 MeV. 

Now we need to specify the viscosity prescriptions used in the hydrodynamical model. We use recent lattice QCD 
calculation results of Meyer[S], for determining C/s. His result indicate the existence a peak of C,/ s near Tc, although 
the height and width of this curve are not well understood. We use the parametrization of Meyer's result given in 
Ref.|22|: 



C 



a exp 



AT 



for T > Tr. 



(25) 



where the parameter a = 0.901 controls the height, AT — Tc/14.5 controls the width of the C/s curve and b = 
0.061. We take the conservative lower bound of the shear viscosity to entropy density ratio rj/s = I/Att in our 
calculations [4 . It is observed that the non-ideal EoS deviates from the ideal case (e — 3P) significantly around the 
critical temperature and around same temperature C/s starts to dominate over rj/s significantly |47j. 
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FIG. 1. Dilepton rate as a function of invariant mass M for different equation of states. Effect of viscosity is not included in 
the hydrodynamical evolution and the distribution functions 



Let us observe the longitudinal pressure of the system given by the equation, 

= P + n - $. (26) 

Since the bulk viscosity contribution 11 < always, 11 and $ together can make Pz negative[22l. (We note that with 
rj/s ^ 1/4tt- as suggested by RHIC experiments, alone is not sufRcient for this condition). The condition = 
defines the onset of cavitation and the time at which it occurs is called cavitation time denoted as Tc- After the 
onset of cavitation the fluid breaks apart and the theory of hydrodynamics is no longer valid to describe the system 
So we can only evolve the hydrodynamics code till in case of occurrence of cavitation instead of r/, where 



V. DILEPTON SPECTRA IN HEAVY ION COLLISION 



The total dilepton spectrum is obtained by convoluting the dilepton rate with the space-time evolution of the heavy 
ion collision. Dilepton rates are temperature dependent and temperature profile is obtained after hydrodynamically 
evolving the system as described in the previous session. In the Bjorken model, the four dimensional volume element 
is given by 

d'^x = cPxTdrisTdr = 7rR\dr]gTdT, (27) 

where Ra — 1.2^1^/'^ is the radius of the nucleus used for the collision (for Au, A — 197). We can calculate different 
differential rates as functions of M, px and y. In this work we will be calculating the rates dN/ {pTdpTdMdy) and 
dN/dMdy; and these dilepton yields are obtained from, 

f dN \ 2 n r""" fl dN \ 

[dM^d^y),,^^^^ - -^-1 ^U^^ 

and are given by 

{prdpXidy) = i^-^^)-^'A r f^^Vs , 

Here tq and ri are the initial and final values of time that we are interested. Generally ti is taken as the time taken 
by the system to reach Tc, i.e.; Tf, but in the case of occurrence of cavitation we must set ti = Tc, the cavitation time, 
in order to avoid erroneous estimation of rates [47j. ynuc is the rapidity of the nuclei used for the experiment. 

Here we note that the dilepton production rates calculated in Section [TTTl correspond to the rest frame of the system. 
So in a longitudinally expanding system, we must replace /o of equation js]) with /g = e"" ^/"^ in equations ( 
With the four momentum of the dilepton parametrised as p" = (raTCOshy,pTCos(j)p,pTsin(j3p,nnTsinhy), where 



(28) 
(29) 



28][29|. 
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FIG. 2. Transverse momentum spectra of dileptons from a viscous QGP calculated at M = 0.525 GeV . The solid line shows 
the dilepton production rate without considering the viscous corrections to the distribution functions. The effect of inclusion 
of viscous corrections due to shear and bulk is shown in separate curves. 



+ M [57] and the four velocity of the medium given by equation ( 19 1 we get, u.p = mTCOsh{y — rjs). Thus using 



the ID boost invariant flow, the factors appearing in the modified rate equation ( 18 1 can be calculated as 



2 2 


4 














T 


r 



(30) 
(31) 



with e = 1/r. 



VI. RESULTS AND DISCUSSIONS 



We calculate the dilepton yields by obtaining the temperature and bulk viscosity as functions of time by solving 
the hydrodynamical equations with relevant initial conditions. For the hydrodynamical evolution of the system we 
use initial conditions relevant for the RHIC experiment, taken from Ref.[15]. Initial time (tq) and temperature (Tq) 
are given as 0.5 fm/c and 310 MeV respectively, whereas ynuc = 5.3. The initial values of viscous terms are taken 
to be zero, i.e.; $(to) — and n(To) — 0. We use the non-ideal EoS (e 



3P ^ 0) obtained from equations (l23|&[24j) 
to close the system. We take critical temperature Tc to be 190 MeV. We will not vary the height or width (controlled 
by the parameter a and AT respectively) of the C,/ s curve in this analysis. These parameters are kept to their base 
values: a — 0.901 and AT — Tc/14.5 throughout this analysis. 

By numerically solving the hydrodynamical equations describing the longitudinal expansion of the plasma ( 20|[22 1, 
we get the temporal evolution profile for T(t), $(t) and n(r). We can evolve the hydrodynamics till the temperature 
of the system reaches critical temperature, i.e.; Tf. 

In order to compare the effect of non-ideal EoS on hydrodynamical evolution and dilepton yields we will compare 
these results with ideal EoS (e — 3P = 0) of a gas of massless quarks and gluons. In this case, within Bjorken flow one 
can solve hydrodynamical evolution equations analytically to obtain the temperature profile as T = Tq (^) s [35]. We 
have studied the temperature profile for both EoSs in our previous work and we found that hydrodynamical evolution 
gets significantly slowed down in case of non-ideal EoS and system spends more time near T^|47j. 

Next we include the another non-ideal effect, viscosity in the calculations. Now we study the longitudinal pressure 
Pz = P + n — $ of the system. It is already seen that in such a scenario, the viscous contribution to the equilibrium 
pressure makes the effective longitudinal pressure of the system zero, triggering cavitation 22 . Hydrodynamics is 
applicable only till Tc in case of occurrence of cavitation instead of Tf. This calculation is presented in our previous 
work[47] in detail and will not be repeated here. We quote the final results here: Eventhough system reaches Tc at 
Tf = 5.5 fm/c only, much before that at Tc — 2.5 fm/c it undergoes cavitation at a temperature 210 MeV. 

Once we get the temperature profile we can calculate the desired dilepton yields as discussed in Section |V] We again 
emphasise that we must be integrating the rates from tq to ti = Tc instead of ri = Tf in the case of cavitation, to 
avoid over-estimation of the yields [l^. From equations ( 28]|29 | we can now calculate the dilepton yields as functions 
of invariant mass M, transverse momentum perpendicular to collision axis px and rapidity y of the dileptons. We 
present all our calculations at the mid rapidity region of the dileptons {y = 0). 
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Fig.[T] shows the dilepton yield dN/dMdy calculated using ideal (massless) and non-ideal EoS. Effects of viscosity 
(both in hydrodynamics and distribution function) are ignored. In this calculation we take Q.bGeV < pt < 2GeV. 
From the figure it is clear that non-ideal EoS yields significantly larger dilepton flux as compared to the ideal EoS. 
At M = 1 GeV, dilepton flux from the non-ideal EoS is about 125% larger than that from the ideal EoS case. This 
behavior can be understood by the fact that system cools slowly in the case of non-ideal EoS allowing a higher 
temperature over a longer period, compared to ideal EoS. It takes almost double the time for non-ideal EoS to reach 
Tc- Since rates are dependent on temperature and an integration over tq and ti — Tf, more dileptons are produced 
in the case of non-ideal EoS. 

It must be noted here that while calculating the particle spectra we use ri = r/ as we have cavitation in the 
system. As we demonstrated in our previous work[3^ particle rates should be integrated upto Tc and if we include 
Tf instead of Tc we will end up having a large over-estimation. In what follows we are presenting the correct particle 
yields by taking into consideration of the effect of cavitation. In Fig. [2] we show the dilepton yield at a fixed low 
invariant mass M — 0.525 GeV as a function of transverse momentum pT of the dileptons. Here we show the effect of 
inclusion of viscous corrections to the distribution function and dilepton rates separately. The solid curve shows the 
case without any viscous corrections to the rates (equation Q): Sf = 0. Now inclusion of shear viscosity corrections 
to the distribution function (denoted by Sf = Sf^j) makes an increase in the dilepton production especially at the 
higher px- This result is in accordance with that of Ref. |51). When we consider only the bulk viscosity corrections 
{6f = Sf(^) we can see that the spectra gets suppressed at the high px regime. This is result is in accordance with 
Refs.[28l [47] where it is shown that effect of bulk viscosity is to suppress particles with high pT- We can see that 
effect of shear and bulk viscosity oppose each other and the total contribution is represented in the curve {Sfjj +(5/^). 
As Pt is increased the corrections are becoming high and the condition fo » Sf may be getting violated. 

Next we consider the same dilepton rate as in previous case, calculated for a high invariant mass M = 1.0 GeV. 
The results are shown in Fig.jsj. The solid curve shows the Sf = case and Sf = Sfr, case is increasing the spectra as 
expected, however corrections to the spectra are becoming increasingly high with high pT- One can see that inclusion 
of modifications due to bulk viscosity alone heavily suppresses the dilepton spectra in this case {Sf = Sfc^). Validity 
of viscous corrections are under scrutiny here as values of both pT and M are high making the corrections diverging. 
Unlike the case for photons, here corrections to the dilepton rates are dependent upon M also (equations (18 30)). 
Our results show that as invariant mass M increases the viscous corrections become large even at smaller px and 
thereby violating the condition fo > Sf. So the validity of the assumption /o > Sf depends on pT as well as M unlike 
the case of thermal photons where only px dependence was there. 

Finally we plot the dilepton rate dN/dMdy as a function of invariant mass of the dileptons in the low pr regime in 
Fig.ji]. We have used PT^i„ = GeV and pt^„^ = 1 GeV in equation (29) while calculating the spectra. It is clear 
that approximation used for finding out viscous modifications in rates are more valid at the regime where both M 
and Pt values are small. 



VII. SUMMARY 



We have calculated the first order corrections arising due to bulk and shear viscosities in dilepton production 
rates from qq annihilation. Viscous corrections to the distribution functions are obtained using the 14-moment 
Grad's method. We have studied the role of non-ideal effects near Tc due to bulk viscosity, EoS and viscosity induced 
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FIG. 4. Invariant mass distribution of mid-rapidity thermal dileptons in RHIC scenario calculated at the low pr regime. Effect 
of viscous corrections to the rate are shown separately. 

cavitation on the thermal dilepton production. Minimal value for shear viscosity rj/ s = l/An is used in in this analysis. 
We have shown that non-ideal EoS taken from lattice results alone can significantly enhance the thermal dilepton 
production. Viscosities also enhance the particle production, however bulk viscosity can induce cavitation much before 
system reaches critical temperature making hydrodynamical description invalid. Thus role of bulk viscosity cannot 
be neglected in particle productions in heavy ion collisions. Finally the viscous corrections in distribution function 
grows in uncontrolled ways if the momentum and/or invariant mass of the dilepton increases and thus the 14- moment 
Grad's method breaks down. 
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